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ABSTRACT 

Temporal solutions are given for the photochemical equa- 
tions describing the distribution of ozone and atomic oxygen 
in an oxygen atmosphere for the equatorial and polar regions. 
Similar calculations are applied to auroral events. These 
calculations indicate a strong dependence of the distribution 
of mesospheric atomic oxygen on the intensity and energy 
spectrum of auroral electrons, which dissociate molecular 
oxygen. It is shown that there can be no significant atmos- 
pheric ozone and atomic oxygen modifications due to soft 
spectrum electron events which appear to characterize most 
bright auroras. The hard spectrum auroral electrons which 
appear in most weak, quiet auroras should cause significant 
increases in the atomic oxygen and ozone concentrations be- 
low 80 km, provided that their flux is more than 10 6 per cm 2 
sec. 
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TEMPORAL VARIATIONS OF MESOSPHERIC OXYGEN AND 
OZONE DURING AURORAL EVENTS 


by 

Kaichi Maeda and A. C. Aikin 
Goddard Space Flight Center 


INTRODUCTION 

The dissociation of molecular oxygen by ultraviolet radiation and charged particles, and the 
subsequent reaction of atomic oxygen with molecular oxygen to form ozone in the earth’s atmos- 
phere was first considered by Chapman (1930). Since then, the chemical equations describing the 
distribution of ozone and atomic oxygen in the mesosphere have been considered by numerous in- 
vestigators including Bates and Nicolet (1950), Dutsch (1961), Barth (1961), Wallace (1962), 
Hesstvedt (1963), Leovy (1964) and Hunt (1964; 1965; 1966 a,b). None of these treatments has con- 
sidered the effect of auroral electrons as a source of molecular dissociation in addition to solar 
radiation. 

Maeda (1962, 1963) has made use of laboratory experiments on the dissociation of oxygen by 
electron impact (Glockler and Wilson, 1932; Massey and Burhop, 1956) to estimate the dissociation 
rate of molecular oxygen due to auroral electrons and protons. Howe/er, more accurate esti- 
mates of the electron impact dissociation cross section as a function of electron energy have been 
reported recently (Lassettre, Silverman, and Krasnow, 1964; Silverman and Lassettre, 1964; and 
Bauer and Bartky, 1965). 

The new data on 0 2 -dissociation by electron impact are employed herein to reevaluate the rate 
coefficient for dissociation by auroral electrons. Also, more rigorous computations are made of 
electron diffusion in the upper atmosphere in which energy loss, coulomb scattering, and the effects 
of straggling are taken into account simultaneously (Maeda, 1965 a,b). The resulting rate coeffi- 
cients are then applied to the calculation of the temporal variation of ozone and atomic oxygen in 
the mesosphere during auroral events. The calculation is accomplished by solving the time- 
dependent photochemical equations for a pure oxygen atmosphere in the 50- to 120-km region. 

The modification of the ionosphere due to the change in concentration of o and o 3 during 
charged particle bombardment is discussed, and the effects of diffusion are presented in Appendix A. 
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DISSOCIATION OF OXYGEN MOLECULES BY ELECTRON IMPACT 


Figures 1 and 2 illustrate the potential energy curves for molecular oxygen and the energy 
level diagram of atomic oxygen, respectively. Excitation of the molecule to repulsive energy levels 
eventually will lead to the dissociation of the molecule, for which the minimum energy is 5.08 ev. 
Photon absorption cross sections in the Schumann -Runge continuum have been measured accurately 
by Watanabe, et al., (1953) and by Ditchburn and Young (1962). Recently, cross sections have been 
determined accurately for dissociation of oxygen by electron impact. 



Dissociative reactions of oxygen 
by electron impact can be categorized 
as follows: 

o 2 + e~ , o + cr , (1) 

-* O + O* + e , (2) 

and 

0 + 0 + + 2e 

- 0~ + 0 + + e > • (3) 

0 + + 0 + + 3eJ 

Reaction 1 is the resonance 
capture of an electron or the so- 
called resonance dissociative attach- 
ment, with the major peak occurring 
at 6.5 ev (Massey, 1950). The total 
cross section for dissociative at- 
tachment, based on recent measure- 
ments by Rapp and Briglia (1964) 
and Rapp, Sharp and Briglia (1964), 
is given in Figure 3. At 6.5 ev the 
maximum cross section for this 
reaction is 1.5 X 10' 18 cm 2 . Reac- 
tion 2, the dissociation of 0 2 into 
two neutral atoms, is of most inter- 
est. The inelastic scattering of 
electrons with incident energies of 
4.5 to 12.5 ev has been studied by 
Schulz and Dowell (1962), and the 
results are presented in Figure 4. 
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The peak occurring around 6.5 ev undoubtedly 
corresponds to Reaction 1. 



Figure 3 — Total cross sections of 0 2 -dissociation by 
electron impacts. Reactions 1, 2, and 3 are indicated 
by c”, ct SL/ and c. , respectively. The previous calcu- 
lation for Reaction 2 (Maeda, 1962, 1963a) is indicated 
by cr M , whereas the calculation for Reaction 2 (Bauer 
and Bartky, 1965) by the classical theory is indicated 
b y ^BB- 


The major peak occurs at 8.4 ev and is due to dissociation from the excited state (b 3 2 u "). 
Dissociation from this level will lead to one atom in the 3 P ground state and one metastable ‘D 
atom (Figures 1 and 2). The peak at 9.0 ev results from dissociation of the *A u state leading to 
two *D atoms. Other peaks up to the first ionization potential at 12.04 ev McGowan et al. (1964), 
are also discernible. Schulz estimated the cross section and found an approximate value of 
1.9 X 10” 18 cm 2 at the 8.4 ev peak. The cross section has also been estimated chemically (Gockler 
and Wilson, 1932; Massey and Burhop, 1956). From an analysis of dc discharges in oxygen gas, 
Thompson (1961) estimated that the cross section for the dissociation of oxygen by electrons of 
energy greater than 14 ev is at least an order of magnitude greater than that for excitation to the 
state. 

B 

The most recent experimental work on the cross section for dissociation of oxygen by electrons 
of energy greater than 20.0 ev has been accomplished by Lassettre (1959), who found that the cross 
section for dissociation by electrons of 500 ev incident energy was 10 percent of the total cross 
section. Lassettre et al. (1964) have shown that the cross section for energy loss at 8.44 ev by a 
500-ev electron is 4 x 10” 16 cm 2 . Most of the dissociations result from this reaction, which is a 
source of O ('d). 


The oscillator strength, derived from the experimental data at 500 ev together with the Born 
approximation was used by Silverman and Lassettre (1964) to derive the dissociation cross section 
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Figure 4— The excitation spectrum of molecular oxygen by electron 
impact between 4 and 12 ev compared with optical absorption bands. 


as a function of incident electron 
energy. This cross section as a 
function of energy is reproduced in 
Figure 3. Also shown is the result 
of a calculation by Bauer and Bartky 
(1965) based on the classical theory 
of Gryzinski (1959) and employing 
the potential energy diagram (Fig- 
ure 1) given by Gilmore (1964). 
Silverman and Lassettre (1964) re- 
gard their calculations as a lower 
bound on the dissociation cross sec- 
tion since only one transition is in- 
cluded. It should be noted that 0 2 - 
dissociation by electron impact is 
far larger than that of photodissocia- 
tion in the same energy range. Fig- 
ure 3 also gives the cross sections 
employed by Maeda (1962, 1963). 

No information is available on 
the dissociation of ozone by elec- 
tron impact. It will be assumed for 
the purposes of calculation that the 
cross section is as large as that 
given for 0 2 . 


Reaction 3 is the dissociative 
ionization of oxygen. Frost and 

McDowell (1958) studied reactions of this type and found that as the colliding electron energy 
increases, not only 0 + , but also O", O, and 0 ++ are produced. 


The cross section, o\, for 0 2 ionization as measured by laboratory experiments is also shown 
in Figure 3 (Tate and Smith, 1932). The dashed curve at the low energy end of these data is an ex- 
trapolation assuming that the ionization cross section vanishes below the 0 2 ionization potential. 


AURORAL DISSOCIATION OF ATMOSPHERIC OXYGEN MOLECULES 
Dissociation by Monoenergetic Electrons 

The energy of most auroral electrons is between 5 and 50 kev, except for occasional hard- 
spectrum auroral electron events (Mcllwain, 1960; Chamberlain, 1961, p. 269; Hultquist, 1964; 
Mann et al., 1963). Since the results of calculations made for monoenergetic electrons can be 
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applied to more general cases of different energy spectra, computations are performed first for 
several cases of monoenergetic electrons. 


Any charged particle entering the earth's atmosphere follows a spiral orbit along a geomag- 
netic line of force. It can be shown, however, that if the magnetic field is uniform and vertical, 
the effect of spiral motion on the diffusion of charged particles in the atmosphere can be disre- 
garded except for the estimation of the horizontal spread in the atmosphere (Maeda and Singer, 
1961; Maeda, 1962; 1963). These charged particles, entering a denser atmosphere, simultaneously 
undergo multiple Coulomb scattering and energy-loss by collisions with air molecules in the up- 
per atmosphere. 

Because of the small mass of electrons, angular spreads of scattered electrons at each colli- 
sion are very large compared with those of protons or heavier particles of the same velocity. 
Energy loss of electrons at each collision fluctuates wildly, causing the so-called straggling effect 
in their residual range, i.e., the penetration depth. Because of the statistical nature of these scat- * 
terings and energy losses, diffusion of electrons into the atmosphere has not been calculated 
rigorously until recently (Spencer, 1959; Rees, 1963; Maeda, 1965 a,b). 

The rate of dissociation of atmospheric oxygen, J d (z, E 0 ), by monoenergetic incident elec- 
trons, E 0 , at the altitude z, can be given by 


J 


de 


r E o 

r E o 

N(z) 

cr d (W) j e (W. E, z)dW 

Jo 

K J 


i(E, E 0 z) dE • dz 


(4) 


where 

N(z) is the number density of molecular oxygen at the altitude z, (cm -3 ); 

o- d (W) is the differential cross section of 0 2 -dissociation by electrons with energy W(cm 2 , 

W in kev) and shown in Figure 3 by °- SL , o- M , and cr BB ; 

j e (W, E, z) is the number of secondary electrons of energy between W and W+ dw, produced by 
electrons of energy E per unit thickness of air at the altitude z (cm -1 , z in km); 

i(E, E 0 , z) is the number of electrons with energies between E and E + dE at the height z, cor- 
responding to the incident energy E 0 , i.e., the differential energy spectrum of 
electrons in energy E and E+ dE at the altitude z, in the atmosphere normalized to 
the incident total flux of monoenergetic electrons of E 0 (cm - 2 sec - 1 ); 

and 

E c (5 ev) is the threshold energy for 0 2 -dissociation. 

The terms inside the square bracket of the integrand (Equation 4) can be written as: 


[ ^dO) j«.( w > E - z ) dW - CT t ' j t ( E ) ' ^( z ) 1 (s) 

•'E 
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where 


cr 


t 



a d (W)dW , 


(6) 


w e (~ 20 ev) is the effective energy for 0 2 -dissociation, 

p( z ) is the density of air (g/ cm 3 ) at the altitude z (km), 

and 

j t (E) is the number of secondary electrons per primary electron with energy E (kev) per 
unit thickness of air (g/cm 2 ). 

This is shown in Figure 5 as a function of E (kev). The curve below 1 kev is derived from the ex- 
perimental data for 0 2 and N 2 (Tate and Smith, 1932); and the residual portion, particularly above 
10 kev, is computed by the following formula: 


j t (E) 


ME) 

V o 


( 7 ) 


where V 0 ^32 ev (Rossi, 1952, p. 47), and k(E) is the rate of ionization loss of electrons with en- 
ergy E in air given by the Bethe formula (Wu, 1960, p. 17). 

The differential spectrum of electrons, i(E, E 0 , z) , at the altitude z in Equation 4 is computed 
from i (e, e 0 , x), the spectrum at the atmospheric depth x, by making use of x-z relation for the 


CIRA-atmosphere (1965). If the ionization loss, 



Figure 5— The number of secondary electrons per unit 
thickness of gas (0 2 and N 2 , in g/cm) per electron of 
energy E (kev). 


Coulomb scattering and the effect of straggling 
are to be taken into consideration simultane- 
ously, then i (e, E 0 , x) can be calculated only by 
the Monte Carlo method. In the present calcula- 
tions, i(E, E 0 , x) is expressed by an analytic 
formula (Appendix B). A comparison of this ex- 
pression with the result obtained by the Monte 
Carlo calculation is shown in Figure 6 for the 
case of vertically incident monoenergetic elec- 
trons; E 0 = 20 kev. The total intensity l(E 0 , <f) 
is also expressed by an empirical formula, (Ap- 
pendix B), which is shown in Figure 7 with re- 
sults obtained by the Monte Carlo calculation 
(for E 0 = 20 kev). In Figure 7, the Monte Carlo 
result of the total energy flux variation with 
penetration depth in air E(f )/E 0 is also plotted 
against € = x/r 0 . 


6 



II 

1 

I 




Figure 7 — Variations of the total intensity of vertical ly 
incident monoenergetic electrons with atmosphere depth. 


r E o 

I(E 0 £) = i (E, E 0 £) 

•'n 


dE 


where £, = X/r 0 , (r 0 = 9.824 • 10 4 g/cm 2 for E Q = 20 
kev). 


Final results of J de (z) cm' 3 sec" 1 are 
shown in Figure 8, as a function of altitude z 
(km) for vertically incident monoenergetic elec- 
trons with an initial energy range of 5 kev £ E 0 £ 1 
Mev. 

By using the approximation (Equation 5 ), the 
right-hand side of Equation 4 is written as: 


Jde( Z ’ E o) dZ = °- t 'N(z)p(z) 


f 

•'E . 


(E) 


• i (E, E q , z) dEdz . (8) 


As can be seen from this expression, J de (z) 
is proportional to o- t , the total cross section of 
0 2 -dissociation by electron impact. The scale 
of the abscissa J de (z) in Figure 8 should be 
read, therefore, by multiplying a factor a, for 
different cases of reactions which are indicated 
by Equations 1, 2, and 3. The factor, a, and the 

effective energy W e (ev) for these reactions and cross sections are shown in Table 1. From this 
table one can see that the application of the Lassettre cross section increases the dissociation rate 


Figure 6— Comparison of i(E, E 0 , x) given by Equation 
B6 (dashed lines) and that of Monte Carlo calculation 
for the vertically incident monoenergetic electrons of 
E 0 = 20 kev (full lines). 
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Figure 8 — Rate of 0 2 -dissociation by monoenergetic 
auroral electrons (with energy E 0 ) as a function of alti- 
tude z (km) for 5 kev < E Q < 1 Mev. 


Table 1 

Multiplication Factors, Cross Sections, and 
Effective Energies for Reactions 1, 2, and 3. 


Reactions 

Cross Section 
(Figure 3) 

at 

W e (ev)* 

1 


0.1 

10 


^SL 

1.0 

80 


°BB 

0.5 

80 

2 

Cr. 

i 

2.0 

150 

3 

cr~ 

<10“ 3 

6 


tMultiplication factors for ( 2 ) as shown in Figure 6. 
tEffective energies for reactions and cross sections. 


by a factor of 10 above that employed by Maeda (1963). The Lassettre cross section value will be 
applied in all subsequent calculations of the dissociation effects of auroral electrons. 


Dissociation by Auroral Electrons with Exponential Spectra 

It has been reported that the energy distribution of the primary auroral electrons, i(E 0 ) > is 
consistent with the following form of exponential spectrum rather than the single power-of-energy 
formula; a monoenergetic distribution, 8 (e-E 0 ), is a fairly good approximation in most cases 
(Mcllwain, 1960; Brown, 1964): 


i( E o) 



(9) 


where E a is a constant (kev). 

The effect of the energy spectrum of the rate of dissociation can be calculated by making use 
of previous results for monoenergetic electrons as follows: 


J D ( z ) 



= y; Yh Jde ( z> E(,i ) exp 

i = 1 



( 10 ) 


where J de (z, E 0 ‘) is given by Equation 8 and shown in Figure 8 as a function of the altitude z (km). 

In Figure 8 E 0 ‘ is taken as a parameter for the range from 5 kev to 1 Mev, and the altitude z is 
indicated in the abscissa. Computation of J D (z) is shown in Figure 9 for E a = 5 and 100 kev. 
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PHOTOCHEMICAL REACTIONS FOR 
MESOSPHERIC OZONE AND ATOMIC OXYGEN 

The amount of ozone above 50 km com- 
prises less than 1 percent of the total ozone 
content in a vertical column of the atmosphere. 

However, a knowledge of its distribution and 
time variation is important in the study of the 
dynamics of the mesosphere, as well as airglow 
phenomena such as the OH-band emission in the 
lower mesosphere (Bates and Nicolet, 1950; 

Wallace, 1962; Hunt, 1966 b). The role of ozone 
in the ionosphere as a possible attaching agent 
to form 0 3 ~ has been discussed in terms of the 
twilight behavior of PCA events (Reid, 1964). 

The altitude distribution and temporal variation 
of atomic oxygen are also important from both the ionospheric and aeronomical standpoints. In 
this section, the time-dependent equations governing the distribution of ozone and atomic oxygen 
will be solved for an oxygen atmosphere. Before investigating the effect of auroral particles, cal- 
culations will be performed for both an equatorial and polar atmosphere under the influence of 
solar radiation. Since the equatorial distribution has been discussed by many investigators, the 
present calculation scheme can be compared with the results of others. 

Equations and Constants 

Principal reactions for ozone formation and destruction are shown in Table 2. The element M 
in the three-body collision is not necessarily oxygen; M is any atom or molecule which balances the 
energies and momentum in the colliding reaction. The number density CM) is, therefore, given by 
the total number density of air. Although the importance of reactions involving hydrogen has been 
treated by Wallace (1962) and shown to be necessary to explain quantitatively the ozone distribution 
more recently derived by Hunt (1966 b), only a pure oxygen atmosphere will be considered in the 
present investigation. 

Corresponding to reactions shown in Table 2, time variations of odd-oxygen allotropes are 
given by the following two equations together with a similar equation for 0 2 : 



Figure 9— Rate of dissociation by auroral electrons with 
energy e~ E o / E a vs altitude z. 


and 


dM 

dt 


d [° 3 ] 

~dt~ - k 


2J 2 [° 2 ] + Js [° 3 ] 


12 [0] [oj [M] - k 13 [0] [Oj - J 3 [0 3 ] 


- k i2 [o] [oj [m] - k 13 [0] [0 3 ] - k n 


2 [0] 2 [M] 


( 11 ) 


( 12 ) 
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Table 2 


Principal Reactions for Atmospheric Ozone 
Formation and Destruction. 

Ozone formation 

O + 0 2 + M — ^>-03 + M 

Ozone destruction 

O + 0 3 kl3 > -20 2 

Atomic oxygen recombination 

O + O + M — + M 

Dissociation of molecular oxygen 

0 2 + h v O + O for \ < 2424A 

0 2 + e — ^>-0 + O + e for E > 5 ev 
Ozone dissociation 

0 3 + hv — ^» 0 + 0 2 for K < 11800A 

0 3 + e ^>-0 + 0 2 + e for E > 1 ev 


where [o], [0 2 ], and [0 3 ] are the number den- 
sities of atomic oxygen, molecular oxygen, and 
ozone, respectively, and [M] is the total number 
density. A summary of the values of the rate 
coefficients k ( , , k 2 2 and k 1 3 employed in the 
present calculations is given in Table 3 together 
with the most recent laboratory determinations. 
The table shows that the reaction rate changes 
with atmospheric temperature, which is a func- 
tion of altitude, latitude, and season. Since the 
main purpose of the present calculation is to de- 
termine the basic character of variation rather 
than to investigate many possible cases, the ef- 
fect of atmospheric temperature will be ignored 
here and discussed separately. 

The total rates of 0 2 -dissociation by solar 
UV- radiations and auroral electrons are given 


Table 3 

Rate Coefficients. 


Reaction 

Present calculation 

Rate 

k n cm 6 sec" 1 
k 12 cm 6 sec" 1 
k 13 cm 3 sec" 1 

3 x 10" 33 
5 x 10" 34 
2.3 x 10" 15 

3 x 10 33 (Reeves et al., 1960) 

8 x 10" 35 exp (445/T) (Benson & Axworthy, 1965) 

5.6 x 10" 11 exp (~285°/T) (Benson & Axworthy, 1965) 


by j 2 and j 3 for 0 2 and 0 3 . The expression for the particle dissociation rate is given by Equation 
8. The photodissociation rate is 


J P i O) “ 2] (™ l ) ’ (is) 

K 

where o-.^ is the dissociation cross section as a function of wavelength for a particular constitutent 
oxygen (i = 2) and ozone (i = 3). The intensity of solar radiation (Q^) outside the earth's atmos- 
phere is represented in units of photons cm" 2 sec" 1 unit wavelength" 1 . The optical depth r ^ repre- 
sents the product of the absorption cross section and total number of absorbing molecules between 
height z and the sun and is a function of time. 

The resulting differential equations are solved with a computer. Since some of the reactions 
are too fast to permit stable integration by the Runge-Kutta method, a more stable, linearized 


10 




finite difference scheme is utilized in the pres- 
ent computation (Appendix C). 

Vertical Distributions of Atmospheric 
Ozone and Atomic Oxygen 

Equatorial Region 

In Figure 10, the diurnal variation of atomic 
oxygen, [O] in cm' 3 , in the equatorial atmos- 
phere at summer solstice is shown for each al- 
titude from 50 to 100 km at 10-km intervals. 
The variation of the vertical distribution of 
atomic oxygen as a function of altitude with time 
as a parameter is given in Figure 11, where in- 
creasing phase and decreasing phase of the 



Figure 11 — Diurnal variation of vertical distribution of atomic oxygen. 



Figure 10 — Diurnal variation of atomic oxygen/cm 3 at 
each level from 50 km with 100 km at 10-km intervals 
in equatorial atmosphere at the summer solstice. 
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Figure 12— Diurnal variation of ozone at each level cor- 
responding to the oxygen-atom variations shown in Fig- 
ure 10, in the equatorial atmosphere at summer solstice. 

90 km used by Hunt (1966) is nearly one order of 
urements (Golomb et al. 1965). 


diurnal variation are presented separately, i.e., 
(a) for morning (increasing phase), (b) for after- 
noon (decreasing phase). Figure 12 shows the 
variations of ozone at each altitude corresponding 
to the oxygen atom variations shown in Figure 10. 
The variation of vertical ozone distribution with 
time as a parameter is given in Figure 13, where 
the increasing phase, which starts after the sun- 
set, and the decreasing phase after the sunrise 
are shown separately in (a) and (b), respectively. 

Comparisons of observations and other cal- 
culations are shown in Figure 14. There is 
reasonable agreement between daytime observa- 
tions and Hunt's calculation for daytime distri- 
bution with respect to the oxygen -hydrogen at- 
mosphere between 50 and 80 km. Although Hunt's 
conclusion, that a pure oxygen atmosphere leads 
to an excess of ozone, can be confirmed by the 
present calculation for the altitude below 80 km, 
his equilibrium (daytime) ozone distribution 
above the 80-km level is too large. It should be 
noted that the oxygen atom concentration above 
magnitude higher than the recent direct meas- 


Qn the other hand, nighttime enhancement of mesospheric ozone, with a maximum of about 70 
km, obtained by present calculation and by Hunt's pure oxygen atmosphere results shows better 




Oj/cm 3 Oj/cm 3 

(a) Increasing phase. (b) Decreasing phase. 

Figure 13— Vertical cross section of temporal variation of ozone shown in Figure 12. 
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Figure 14 — Comparison of observed values and theoretical results of ionospheric ozone distributions 
at the daytime equilibrium and nighttime enhancement. 


agreement with observations that oxygen-hydrogen atmosphere results. Since the data of Reed and 
Scolnik require some corrections for the spectrum-line widths with altitude,* detailed comparison 
with data seems premature. It should be noted, however, that Mikirov’s observation indicates far 
larger nighttime enhancement than present estimation and rather agrees withNicolet's and Horiuchi's 
estimations. Further nighttime mesospheric ozone measurements are necessary to discuss these 
differences. It has been noted that below 80 km, the nighttime ozone concentrations are unaffected 
by hydrogen reactions, since in the absence of sunlight the free hydrogen concentration falls to in- 
significant values. Therefore, unless there is significant dissociation of water vapor by electron 
impact, a pure oxygen atmosphere should be sufficiently good approximation for use in describing 
auroral events occurring at night. 


*Edith Reed, Private communication, 1966. 
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From Figures 10, 11, 12, and 13, the re- 
laxation time of the atmospheric oxygen atom, 

Tj (sec), and that of ozone r 3 (sec) can be ob- 
tained for the temperature-independent rate 
coefficients of Table 3. The time constants are 
compared in Figure 15, with Hunt's results 
(1964). 

Polar Region 

The distribution of constituents in the polar 
winter mesosphere is particularly difficult to 
describe because of prolonged twilight condi- 
tions. In Figure 16(a) and (b), corresponding to 
the solar declination, 8=0 degrees (equinox) 
and 8 = -10 degrees, respectively. On this basis, one would expect very low atomic oxygen con- 
centration in the winter polar mesosphere, if there is no subsidence and horizontal transport from 
the sun-lit latitudes. Kellogg and Schilling (1951) considered adiabatic heating to be due to air- 
subsidence at least at the rate of 1 km per day in an effort to explain the increase in temperature 
of the winter mesosphere, which is warmer than that in summer. Later Kellogg (1961) suggested 
that atomic oxygen created at heights above 80 km would be carried polewards to the dark hemis- 
phere, where subsidence of the oxygen would increase the heating rate by the exothermic (5.08 ev) 
three-body recombination reaction. Young and Epstein (1962) have shown that rates of subsidence 
of 0.2 cm per sec are sufficient to account for the observed heating provided that the atomic oxygen 
concentration at 115 km is 1.5 x 10 11 cm -3 or greater. On the basis of this density, it was shown 
that atomic oxygen densities of the order of 10 12 cm' 3 can be obtained at 70 km. The value of the 
rate coefficient which was chosen for k 13 in the calculation is, however, one order of magnitude 
larger than that utilized by Hunt (1966). It is also possible for diffusion to be active in redistri- 
buting the atomic oxygen concentration and ozone formed in the upper mesosphere and in the lower 
thermosphere; the effect of diffusion is estimated in Appendix A. 



Figure 15 — Time constants of atmospheric oxygen-atom 
r 1 (sec) and of ozone, t 3 (sec), and comparison with 
Hunt's (1964) results. 


MESOSPHERIC OXYGEN DISTRIBUTION DURING AURORAL EVENTS 

In Figure 9, the dissociation rates per incident electron have been given for two expo- 
nential spectra for e-folding energies, E a = 5 and 100 kev. These spectra represent the 
limits for soft and hard spectrum of observable auroral electrons and can be used to establish 
the particle flux necessary to modify the atomic oxygen distribution in each instance. Such 
an estimate can be made by employing the time-independent solution of Equation 12 and 
solving for the flux necessary to overcome the decay terms. As can be seen from Figure 9, a 
spectrum with E a = 5 kev can have no effect at altitudes of 80 km and below. However, if a flux 
is more than 5 x 10 s electrons per cm 2 sec, the production rate can overcome recombination effects 
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Figure 16(a)— Temporal variation of ozone and atomic oxygen in the polar hemisphere at the geographic 
latitude \ = 70 degrees and solar declination S = 0 degrees (equinox). 


at 70 km for the E a = 100 kev case. It is, therefore, the hard spectrum events that are of impor- 
tance in modifying the atmosphere below 80 km, where atomic oxygen has a short lifetime (Figure 
15). Most bright auroras, however, are not of this hard spectrum type, but are characterized by 
the soft spectrum with E a ’s of 5 to 10 kev and by the fluxes of 10 9 to 10 12 electrons per cm 2 sec 
(Mcllwain, 1960; Brown, 1966). Such events are not important unless the duration is of the order 
of several hours, which might affect the atomic oxygen distribution above 80-km level, since the 
lifetime of atomic oxygen is long above these altitudes. Hard spectrum events of the type sufficient 
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Figure 16(b) — Temporal variation of ozone and atomic oxygen in the polar hemisphere at the geographic 
latitude A- = 70 degrees and solar declination S = -10 degrees. 

to affect the atomic oxygen concentration have been reported by Bailey, et al. (1966) on the basis of 
forward -scatter radio-propagation data and balloon measurements of auroral Bremsstrahlung (Brown, 
1966). Such events, referred to as relativistic electron precipitation events, or REP's appear to be 
characterized by E a 's between 60 and 150 kev and fluxes between 10 6 and 2 X 10 4 electrons per cm 2 
sec. There is also satellite evidence of hard spectra auroral electrons (Mann, et al., 1963). 

In the present calculation, Equations 11, 12, and 13 are solved by the method described previ- 
ously (and shown in Appendix C) for a latitude of 70 degrees during equinox. In addition to the small 
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solar photodissociation rate, particle fluxes are superimposed. As stated previously the cross 
section of 0 3 dissociation by electron impact is not known. However, the dissociation cross sec- 
tion of ozone (a triatomic molecule) by electron impact should be larger than that of the diatomic 
molecule. Therefore, two values have been assigned to the ozone dissociation rate; i.e., J 3 = 10 2 J 2 
and J 3 = J 2 . It will be shown that the effect of auroral events on the temporal and spatial varia- 
tions of atmospheric ozone and oxygen are not drastically affected by the assumption of the 0 3 - 
dissociation rate unless it is assumed to be unreasonably small. 

The following eight cases have been considered in the calculations: flux (intensity) of imping- 
ing auroral electrons i 0 = 10 11 to i 0 = 10 s (cm -2 sec' 1 ); energy spectrum of auroral electrons 

i ( E o) = i o ex P (~ e ^) ’ ( 14 ) 


where E a = 100 kev (hard spectrum) and E a = 5 kev (soft spectrum); and the ratio of dissociation 
rate of 0 2 and of 0 3 , J 3 = J 2 x 10 2 and J 3 = J 2 . Duration of auroral bombardment was assumed 
to be 10 minutes and initiated at midnight for all calculations. 


It was found that the soft spectrum case had essentially no effect even for large fluxes of in- 
cident electrons. Hard spectrum events of even short duration (1 minute or so) can give rise to 
significant enhancements at low altitudes. For instance, for a flux of 10 5 electrons per cm 2 sec 


a density 10 5 cm' 3 atomic oxygen can be main- 
tained at night at 60 km and for the same flux 
the density at 70 km, can reach densities of 10 8 
cm -3 . The resulting atomic oxygen and ozone 
distributions for the hard spectrum case (E a 
= 100 kev) are shown in Figures 17 (a through 
c) and 18 (a through c) for J 3 = J 2 X 10 2 and 
J 3 = J , respectively, for the extreme case of 
i 0 = 10 11 cm' 2 sec' 1 . 

From these figures, it can be seen that the 
atomic oxygen concentration is altered signifi- 
cantly below 80 km. In fact, there is an in- 
crease of concentration of five orders of mag- 
nitude at 70 km, which decays with time constant 
Tj indicated in Figure 15. The increase in 0 3 
is most significant around the 60 to 70-km level 
and is almost independent of the assumed ozone 
dissociation rate unless this is chosen to be 
negligibly small. From the foregoing, it is 
evident that for a hard spectrum auroral elec- 
tron (E a = 100 kev) and large incident flux 



TIME ( hours) 

Figure 17(a) — Variation of the O and 0 3 concentration 
during and after the hard spectrum electron precipita- 
tion at 60 and 70 km, for J 3 = J 2 xl0 2 andi 0 = 10" 
cm ' 3 sec -1 . 
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Figure 17(b)— Variation of the O and 0 3 concentration 
during and after the hard spectrum electron precipita- 
tion at 80 and 90 km for J 3 = J 2 x 10 2 and l Q = 10 11 
cm -3 sec -1 . 
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Figure 18(a)— Variation of the O and 0 3 concentration 
during and after the hard spectrum electron precipita- 
tion at 60 and 70 km for J 3 = J 2 . 
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Figure 17(c)— Variation of the O and 0 3 concentration 
during and after the hard spectrum electron precipita- 
tion at 100, 110, and 120 km for J 3 - J 2 x 10 2 and i 0 
= lOH cm -3 sec"’ . 



Figure 1 8(b)— Variation of the O and 0 3 concentration 
during and after the hard spectrum electron precipita- 
tion at 80 and 90 km for J, = J.,. 






(i 0 = 10 11 cm' 2 sec" 1 ), the profiles of atomic 
oxygen and, to a lesser extent, of ozone are 
substantially modified during and after the 
period of electron precipitation. Since the dif- 
fusion time of dissociated oxygen is long at 
these altitudes, the enhancement is confined 
within the horizontally narrow domain of par- 
ticle precipitation (within bright auroras). 

From the present calculation, it is seen 
that the effect of O a dissociation is significant 


Ki- 


lo' 


10 7 


10 c 


F [O] 100 km j 


r [O 3 ] 100 km 




10 4 


10 J 


b- [o] 110 km 4 


f [O 3 ]ll0km | 


_L 


E [O] 120 km j 


\ [O 3 ] 120 km 


12 


14 


16 12 14 16 

TIME (hours) 


12 


14 


16 


Figure 18(c)— Variation of the O and 0 3 concentration 
during and after the hard spectrum electron precipita- 
tion at 100, 1 10, and 120 km for J 3 = J 2 . 


only below the 80-km level for both processes 
of dissociations, i.e. by auroral particles as 
well as by solar UV-radiation. Since no auroral 
electrons with energy less than 100 kev pene- 
trate below the 80-km level (Figure 8), the rate 
of auroral dissociation depends strongly on the 
amount of high energy electrons or on the spec- 
trum of auroral electrons. In other words, if 

the spectrum is soft as seen in most bright auroras (Mcllwain 1960, Brown 1966), no effects of 
auroral dissociation are detected. On the other hand, if the spectrum is hard, an increase of 
atomic oxygen below the 80-km level, together with a subsequent increase of mesospheric ozone, 
can occur. It should be noted, however, that hard spectra of auroral electrons appear generally in 
weak displays (Brown 1966). Durations should be, therefore, more than hours before the effects of 
auroral dissociation of atmospheric oxygens can be observable. 


If local variations of auroral activities are considered, the intensity of the auroral electron, 
and its spectra change drastically within milliseconds and several hundreds meters of horizontal 
scale (O'Brien, 1964; Brown, 1966; Mozer and Bruston, 1966). The time constants of oxygen al- 
lotrope variations above stratosphere are, however, at least in the order of minutes (Figure 15). 
The intensities and spectra of auroral electrons to be compared to the present calculation are, 
therefore, those averaged over several minutes and many hundreds of kilometers in horizontal 
scale. In this respect, auroral effects on atmospheric oxygen are quite similar to those deduced 
from the auroral absorption measurements by riometer as far as time and horizontal scales are 
concerned. 


Finally, the following can be regarded as one of the important differences between the effects 
of solar UV-radiations and those of auroral particles on the atmospheric oxygen; i.e., in the case 
of solar UV-dissociation, the increase of atomic oxygen is always accompanied by the decrease of 
ozone, while both oxygen atoms and ozone increase in the case of auroral dissociation. This is 
essentially due to the steep spectrum of solar UV-radiation in which the intensity of the Schumann- 
Runge region is nearly three orders of magnitude less than those ozone -destroying Hartley 
regions. 
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According to Akasofu (1964), the input energy flux of auroral particle is of the order of 80 ergs 
per cm 2 sec (the total energy deposit-rate, 8 X 10 16 ergs per sec, divided by the total cross- 
sectional area exposed to incoming auroral electrons, 10 15 cm 2 ). This figure is comparable to the 

O 

energyflux of solar UV-radiation with wavelength shorter than 2000A (Hinteregger, 1965). It should 
be noted, however that the energyflux of bright aurora exceeds this figure more than several orders 
of magnitude, fluctuating within the order of milliseconds. 


IONOSPHERIC EFFECTS OF OZONE AND ATOMIC OXYGEN ENHANCEMENT 

When an auroral electron interacts with the atmosphere, it loses most of its energy through 
excitation and ionization. The resulting ionization leads to increased radio wave absorption dur- 
ing auroral electron precipitation, and in the case of proton flares, polar cap absorption events. 

An adequate description of these events depends on a knowledge of the electron distribution, which 
is dependent on the negative ion distribution in the D region. The species of negative ion is un- 
known at present. However, it is usual to assume 0 2 ~ although O" and N0 2 ~ have been suggested. 
The negative ion is formed initially by electron attachment. The species of negative ion may then 
be changed by charge transfer reactions, or the ion may be destroyed by a variety of processes. 
These include the very important associative detachment reaction 


0 2 " + O - 0 3 + e . (15) 

This reaction has figured in practically every theoretical treatment of the lower ionosphere 
both for the normal diurnal variation Aikin (1962) and for disturbed events such as PCA's (Whit- 
ten and Poppoff, 1962; Reid, 1966). Associative detachment was also utilized to explain auroral 
absorption events by Aikin and Maier (1963) who suggested that a modification of the atomic oxygen 
distribution by the incoming energies particle might change the negative ion density profile during 
such events. 

Fehsenfeld et al. (1966) have recently measured the rate of Equation 15 and obtained a value 
of 3 x 10" 10 cm -3 sec -1 which is more than sufficient to make this process dominant above photo- 
detachment during the day. If only 0 2 " is considered, the negative ion distribution will be controlled 
entirely by the atomic oxygen concentration. However, the possibility that 0 2 " and O” may charge- 
exchange with neutral species to form ions such as 0 3 " and N0 2 " must be given consideration. 
Ferguson* has measured the rate of the reaction 

° 3 " + °°2 ’ C 0 3 ~ + °2 ’ 

and found it to have a rate of 3 x 10" 10 cm -3 sec" 1 so that ions of this species are a possibility as 
well as 


0 3 " + NO - NO~ + 0 . 


*E. E. Ferguson, Private communication, 1966. 
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It is evident that the distribution of O and 0 3 can greatly influence the ionosphere wherever 
negative ions are of importance. The possibility exists that the negative ion distribution and per- 
haps the species is quite different between the normal D region, that affected by hard spectrum 
auroral electrons and PCA protons, simply because the distribution of ozone and atomic oxygen 
differs in the lower mesosphere during these circumstances. 


CONCLUSIONS 

The time -dependent photochemical equations describing the vertical distribution of ozone and 
atomic oxygen in a pure oxygen atmosphere have been solved for the equatorial and polar regions. 
The extended periods of twilight play an important role in the polar regions. Continuous sunlight 
can lead to an excess of atomic oxygen at altitudes where the atomic oxygen lifetime is greater 
than one day. For winter periods when the insolation is absent, appreciable atomic oxygen con- 
centration cannot be maintained at altitudes below 80 km without a downward flow mechanism, such 
as subsidence or diffusion, or without a production mechanism, such as a large flux of energetic 
electrons. 

It has been shown that the cross section for molecular oxygen dissociation by low energy elec- 
trons is quite large (Figure 3). Thus secondary electrons created during auroral particle precipi- 
tation into the atmosphere can modify the atomic oxygen and ozone distribution. If the energy spec- 
trum of the auroral primary electrons is very hard such that the e-folding energy E a is 50 kev or 
greater, mesospheric ozone and atomic oxygen concentrations around 60 to 70 km increase signifi- 
cantly (Figures 17 through 20). The majority of events are not of this type but are generally very 



(a) During electron precipitation 
(increasing phase). 


(b) After the cessation precipitation of 
J 3 = J 3 x 10 2 and i 0 = 10 11 cm" 2 sec" 1 . 


Figure 19 — Variation of the vertical distribution of atomic oxygen due to the hard spectrum electron precipitation. 
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Figure 20— Variation of the vertical distribution o 
ionospheric ozone during and after hard spectrum elec 
tron precipitation for J 3 = J 2 x 10 2 and i Q = 10 1 
cm -2 sec -1 . 

emissions lead to 0 2 -ionization rather than ( 
not affect present results. 


soft and rather monoenergetic with an energy 
range of 10 to 50 kev. For events of this type, 
there is no modification of the atmosphere un- 
less the event is of such a long duration as to 
increase the atomic oxygen concentration above 
90 km. Therefore, significant enhancement of 
mesospheric ozone during auroral activity 
(Murcray 1957) seems very unlikely in general, 
although it may occur for special events. Most 
of the UV emissions from excited nitrogen 
molecules originating from auroral electron 
i bombardments are in the far ultraviolet re- 
gions which exceed the ionization threshold of 
molecular oxygen. Therefore, the auroral UV 
2 -dissociation (Green and Barth 1965) and will 
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Appendix A 


Effects of Diffusion 


The diffusion of neutral atoms and molecules in the upper atmosphere can be described in 
terms of molecular diffusion and eddy diffusion, both of which depend on the vertical gradients of 
densities, pressures, and temperatures in the atmosphere. Calculations on these effects have been 
made by several investigators (Mange 1954; Colegrove et al., 1965; Shimazaki, 1966). Numerical 
computations which include these effects are not impossible, but they create additional complica- 
tions and increase machine time. Before including such processes, the effects of diffusion can be 
estimated. Since the products of 0 2 -dissociation have a finite lifetime in the atmosphere, which 
changes with altitude, one can define the effective diffusion length of these elements by 

( z ) = [ D i ( z ) ' T i ( z )] 1/2 - 


where 

D. ( z >'s are the diffusion coefficients for 0 ; at the altitude z km, in cm 2 sec -1 , and 

r. ( z )'s are the lifetime of Oj in the atmosphere at the altitude z km, in second, which are 
shown in Figure 15. 

It is known that above a certain altitude (turbopause), the molecular diffusion coefficient ex- 
ceeds the eddy diffusion coefficient and increases further with altitude. Since the present calcu- 
lation is an order of magnitude estimation, (z)'s are taken from the corrected value given by 
Johnson and Wilkinson (1965). If (z) is less than the local vertical scale for gradient hj (z) of 
the element 0 ; at z in km, the effects of diffusion can be regarded as unimportant; while if -t. (z) 
exceeds hj (z), diffusion is important to consider in finding the vertical distribution of elements 
Oj in the atmosphere around the altitude z. The results are shown in Table Al. From this table, 
we can conclude that the diffusion is not important to determine the vertical distributions of atomic 
oxygen and ozone below 90 km but cannot be disregarded above the 100-km level. As can be seen 
from Figure 15, the effect of diffusion is not as important for the ozone distribution as for that of 
atomic oxygen. 
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Table A1 


The Effective Diffusion Length, / t i (km) and the Local Vertical 
Scale for Gradient, h i (in km) in the Mesosphere, Where 
Suffix i = 1, 3 stand for Atomic Oxygen and Ozone, Respectively. 


Altitude 
z in km 

(oxygen atom) 

hi 

(oxygen atom) 

^3 

(ozone) 

h 3 

(ozone) 

50 

0.06 

1 

0.22 

5.0 

60 

0.1 

3 

0.50 

6.5 

70 

0.35 

6 

0.27 

5.0 

80 

1.8 

12 

0.18 

4.5 

90 

8.7 

21 

0.08 

4.0 

100 

120 

CO 

0.054 

3.5 

110 

3600 

13 

0.036 

3.0 

120 

>10 5 

20 

0.020 

3.0 



Appendix B 


Analytic Formula for i(E, Eg, x) 


The number of electrons whose energy is between E and E + dE (in kev) at the atmospheric 
depth x (g/ cm 2 ), corresponding to monoenergetic vertically incident electrons with energy 
E 0 , i (e, E 0 , x) has been calculated by means of the Monte Carlo method (Maeda, 1965). This 
can be approximated by 


Ae ~ a (^o) b 

1 (^’ ’ x ) “ N (E 0 , x) * o ’ x ) ’ 


(Bl) 


where 


and 


and 



(B2) 

(B3) 

(B4) 
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The empirical formula for a range r 0 (g/cm 2 ) of electrons with initial energy E 0 (kev) in air 
shown by Equation B3 (Grun, 1957) and E m given by Equation B4 corresponds to the residual energy 
of electrons at the depth x in air. 

Since the integrand i(y, 1, x) is independent of E 0 , the normalization factor of Equation Bl, 
n(e 0 , x), is written as: 


n(e 0 . x) 


e„ 



(y, x) dy 


(B5) 


where 


i o (y- x) 


f d-y) 1 

"GO 

r l 1-y y‘ f ’l 

k-yX'i 

exp 



(B6) 


Numerical values of Equation B6 for E 0 = 20 kev are shown in Figure 6, with comparison to 
Monte Carlo results (Maeda, 1965 a,b). 

Another term in Equation Bl, Ae ^ + ^°' ) , is an empirical expression of total intensity of verti- 
cally incident monoenergetic electron E 0 , at the depth x = • r 0 , which also has been computed by 

the Monte Carlo method (Maeda, 1965). Comparisons of numerical values of this factor and Monte 
Carlo results are shown in Figure 7. 
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Appendix C 


Integration of the Differential Equations by Finite Difference Method 


The differential equation (Equation 11) can be replaced with the following finite difference 
scheme; 


0 3 = ° 3 At — = 0,<"> [k la -M-0 2 <" +,) -k 13 •0 3 <n+1 >] - J 3 -0 3 <" +1 > , (Cl) 

where o it cr stand for [o.j and d [o ; ] j dt in Equation 11, and the superscripts n, n + 1 refer to the 
time t n and t n+1 = t n +At , respectively. 

Similarly, Equation 12 can be replaced by 


0 < n+ 1 ) - o 

Oj - - 1 At — = 2J 2 ' 0 2 < n + l > + J 3 • 0 3 < n+1 > + Oj< n > [- 2k u ■ M • 0/ n+1 > 


-k 12 -M-0 2 <" +1 )-k 13 -0 3 <" + 1 >] . (C2) 


Corresponding to the reactions listed in Table 2, time variation of 0 2 also can be written as; 


o 2 < n+1 > -0 2 < n > 
At 


= 0/"’ jk u M •0 1 (n + 1 > + 2k 13 0 3 < n+1) ] - 0 2 < n + 1 > • [kjjM •0 1 (n + 1 > + J 2 ] + J 3 • 0 3 < n+1) .(C3) 


This equation can be also derived by Equations 11 and 12 with a condition for conservation of oxygen 
allotrope at each layer; i.e., 


3 [0 3 ] + 2 [0 2 ] + [Oj] = const. 


or 


[0» +30 3 ] 


(C4) 
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It should be noted that Equations Cl, C2, and C3 are written in linear form which can be simpli- 
fied by writing in matrix, i.e., 


/o/ n) \ 

° 2 <n) 

\° 3 (n 7 


( 1 +Oj (n) • 2k n MAt 
-0/ n) kjj MAt 
° 


(0,< n > k 12 M - 2J 2 ) At 
1 + (0Wk 12 M+J 2 ) At 
-0/ n) k 12 MAt 


lO/ n)k i3 -J 3 ) At 
- (0,< n> 2k 13 + J 3 ) At 

1 + ((>,<"> k l3 + J 3 ) At 


/OjC n+1 ) 

0 2 (" + i) 

\ 0 3 (n+1) 


\ 

/ 


By means of matrix inversion O i < n+1) 's are easily obtained by knowing Oj (n) 's, starting from 
certain initial conditions. As can be seen from Equations Cl and C3, replacement of original dif- 
ferential Equations 11 and 12 by the finite difference schemes is not unique. It should be noted, 
however, that partial replacements such as o i (n+I) • o. (n + 1) by Oj (n) • 0. (n+1) are essential to retain 
the linearity of the scheme by which integration stability in numerical integrations is to be 
achieved easily. 

The present method works for At = 20 sec for altitudes of 50 to 150 km and seems progres- 
sively unstable below 50 km for At = 20 sec, which can be extended, however, by reducing At. 


34 


NASA-Langley, 1967 


13 



"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

—National Aeronautics and Space Act of 1958 
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